Statistics of reflection eigenvalues in chaotic cavities with non-ideal leads 
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The scattering matrix approach is employed to determine a joint probability density function of 
reflection eigenvalues for chaotic cavities coupled to the outside world through both ballistic and 
tunnel point contacts. Derived under assumption of broken time-reversal symmetry, this result 
is further utilised to (i) calculate the density and correlation functions of reflection eigenvalues, 
and (ii) analyse fluctuations properties of the Landauer conductance for the illustrative example of 
asymmetric chaotic cavity. Further extensions of the theory are pinpointed. 
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Introduction. — At low temperatures and voltages, a 
phase coherent charge transfer through quantum chaotic 
cavities is known to exhibit a high degree of statistical 
universality^. Even though the transport through an in- 
dividual chaotic structure is highly sensitive to its micro- 
scopic parameters, the universal statistical laws emerge 
upon appropriate ensemble or energy averaging proce- 
dure. The latter efficiently washes out all system-specific 
features provided a charge carrier has stayed in a cavity 
long enough to experience diffraction 3,4 . Quantitatively, 
this requires the average electron dwell time td to be in 
excess of the Ehrenfest time te that defines the time scale 
where quantum effects set in. 

In the extreme limit td 3> te , the statistics of charge 
transfer is shaped by the underlying symmetries^ of a 
scattering system (such as the absence or presence of 
time-reversal, spin- rotational, and/or particle-hole sym- 
metries). For this reason, a stochastic approach^ based 
on the random matrix theory^ (RMT) description! of 
electron dynamics in a cavity is naturally expected to 
constitute an efficient framework for nonperturbative 
studies of the universal transport regime. Indeed, a stun- 
ning progress was achieved in the RMT applications to 
the transport problems over the last two decades. Yet, in- 
tensive research in the field^ left unanswered many basic- 
level questions. One of them, regarding the statistics 
of transmission/reflection eigenvalues in chaotic cavities 
coupled to the leads through the point contacts with tun- 
nel barriers, will be a focus of this Letter. 

Supported by the supersymmcry field theoretic 
technique 9 as well as by recent semiclassical studies^, 
the RMT approach to quantum transport starts with the 
Heidelberg formula for the scattering matri>*i^ 

S(e F ) =1 N - 2mtW + (£ F t M -0C + iTrWW^^W (1) 

of the total system comprised by the cavity and the leads. 
Here, anMxM random matrix 5C (of proper symmetry, 
M — > oo) models a single electron Hamiltonian whilst 
an M x N deterministic matrix "W describes the cou- 
pling of electron states with the Fermi energy e F in the 
cavity to those in the leads; N = + n K is the to- 
tal number of propagating modes (channels) in the left 
(til) and right (na) leads. Equation (fTJ) refers to chaotic 
cavities with sufficiently large capacitance (small charg- 



ing energy) when the electron-electron interaction can be 
disregarded-!.. Throughout the paper, only such cavities 
are considered. 

Landauer's insight 11 that electronic conduction in 
solids can be thought of as a scattering problem makes 
the N x N scattering matrix S(sf) a central player in 
statistical analysis of various transport observables. In 
the physically motivated M — > oo scaling limit, its dis- 
tribution, dictated solely by the symmetries of the ran- 
dom matrix IK, is well studied for both normal 1 — and 
normal-superconducting 13 chaotic systems. In the for- 
mer case, the distribution of S(e_f) is described by the 
Poisson kerne l 12 ' 14 



P /3 (S) oc [det(l A r-S§ t )det(l A r-SS t )] 



/3/2-1-/3AT/2 
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Here, /3 is the Dyson indexpi accommodating system sym- 
metries: §(sf) is unitary symmetric for f3 = 1, unitary 
for (3 = 2, and unitary self-dual for (3 = 4. All relevant 
microscopic details of the scattering system are encoded 
into a single average scattering matrix 



S = (MAIjv - n 2 W f W) (MAI 



JV 
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where A denotes the mean level spacing at the Fermi 
level in the limit M — > oo. The TV eigenvalues 7 — 
diag({-\/l — Tj}) of S characterise^ couplings between 
the cavity and the leads in terms of tunnel probabili- 
ties Tj of the j-th electron mode in the leads. The cele- 
brated result Eq. ([2]) , that can be viewed as a generalisa- 
tion of the three Dyson circular ensembles^, was alterna- 
tively derived through a phenomenological information- 
theoretic approach reviewed in RefJ^. 

Unfortunately, statistical information accommodated 
in the Poisson kernel is too detailed to make a nonpertur- 
bative description of transport observables operational. 
It turns out, however, that in case of conserving charge 
transfer through normal chaotic structures, it is suffice 
to know a probability measure associated with a set T 
of non-zero transmission eigenvalues {Tj £ (0, 1)}; these 
are the eigenvalues of the Wishart-type matrix tv , where 
t is the transmission sub-block of the scattering matrix 



t' 



n R xn L 



n R xn R 
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Owing to this observation, the joint probability density 
function Pp(T) emerges as the object of primary interest 
in the RMT theories of quantum transport. 

Surprisingly, our knowledge of the probability measure 
Pg(T) induced by the Poisson kernel [Eq. ©] is very lim- 
ited, being restricted to chaotic cavities coupled to exter- 
nal reservoirs via ballistic point contacts^ 6 ("ideal leads" ). 
In this, mathematically simplest case, the unity tunnel 
probabilities Tj = 1 make the average scattering ma- 
trix S vanish, giving rise to the uniformly distributed 1 ' 
scattering matrices which otherwise maintain a proper 
symmetry^. In the RMT language, this implies that scat- 
tering matrices belong to one of the three Dyson circular 
ensembles 1 ^. 

As was first shown by Baranger and Mello 1 ^, and by 
Jalabert, Pichard, and Beenakker— , the uniformity of 
scattering matrix distribution induces a nontrivial joint 
probability density function of transmission eigenvalues 
{Tj} of the form 2 ! 



p (/3) (t)c< |A£(T)inr. 



P/2-1+0V/2 



(5) 



min(nL, iir) is the number of non-zero eigenvalues of the 
matrix tt\ whilst A„(T) is the Vandermonde determi- 
nant A„(T) = Ylj <k (Tk-Tj). Equation (5) is one of the 
cornerstones of the RMT approach to quantum trans- 
port. 

From ballistic to tunnel point contacts. — The restricted 
validity of Eq. ([5]), that holds true for chaotic cavities 
with ideal leads, is hardly tolerable both theoretically 
(an important piece of the transport theory is missing 2 ^) 
and experimentally (chaotic structures with adjustable 
point contacts, including tunable tunnel barriers, can 
by now be fabricated 2 ^). In this Letter, a first system- 
atic foray is made into a largely unexplored territory 
of non-ideal couplings. In doing so, we choose (for the 
sake of simplicity) to lift a point-contact-ballisticity only 
for the left lead which is assumed to support prop- 
agating modes characterised by a set of tunnel prob- 
abilities Tl = (Tx, . . . ,r„ L ); the right lead, support- 
ing riR. > til open channels 24 , is kept ideal so that 
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n L + l) ■ 



1, 



Assuming that the 



Here, v = \n\, — Ur| is the asymmetry parameter, n 



time-reversal symmetry is broken (/3 = 2), we shall show 
that the joint probability density function ) (-R) of 
reflection eigenvalues {Rj = 1 — Tj} equals 2 ^ 



- P (7L|0)(-Rl I 



! ) 



det JV (l„ L -^ ) 



A nL (il)det (j!fc)e(1 , nL) [ 2J Fi(n R + l,n R + l;l;7? J R fe )l J](l-i? i ) I/ . (6) 
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Here, 7^ = l nL — lY is a set of coupling parame- 
ters characterising non-ideality of the left lead in terms 
of associated tunnel probabilities, N = til + «r, is the 
total number of open channels in both leads, c„ L! „ R is 
the inverse normalisation constant, 



(n L + n R )l 



K0 



n 2 



n L \n K 



(n R + j)!(n R - j)! 



(7) 



whilst p F q is the Gauss hypergeometric function. 
The (biorthogonal 27 ) ensemble of reflection eigenvalues 
Eq. (IB]) is our first main result 26 . Before outlining its 
derivation, let us discuss the implications of Eq. © for 
a nonperturbative statistical description of both spectral 
and transport observables in quantum chaotic cavities. 
Statistics of reflection eigenvalues. — The first immedi- 



ate consequence of Eq. ([6]) is the determinant structure of 
the p-point correlation function of reflection eigenvalues: 

P(n L ,n R )(Rl, ■ ■ ■ ,Rp) = det(j,fe) e (i, p ) [if(„ Li „ R )(i?j,i4)] . 

.(8) 

Defined in a standard manner—, it is entirely determined 
by the two-point scalar kernel -f^(n L ,n R ) {R, R'), that can 
straightforwardly be calculated 2 ^ by applying the ideas 
exposed in Refill. In terms of the "moment function" 



M, 



(9) 



where \ 1 2 ) = k~ x 3 F 2 (n R + l,n R + 1, k; 1, £;+ 1; 7 2 ), 
the scalar kernel is given by a finite sum: 



Kin^nioiR, R') = «l!c„ Li „ r 



det 



( f" 2 \ ^ [{l-R){l-R'W /2 



An L (7E) 



(10) 



x ^^^nR + l^R + Ll^^det^IM^^ 2 )]^,...^!; (tf)*" 1 ; [Mf( 7 ])] H+1 , 
3=1 
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FIG. 1: Probability density function /(i,n R )(<?; T) for Lan- 
dauer conductance plotted for tir = 5 and various tunnel 
probabilities: V = 0.99 (solid line), F = 0.8 (long-dashed line), 
r = 0.6 (dashed line), F = 0.4 (dotted line), and T = 0.2 (dot- 
dashed line). For the ideal point contact (r = 1), the curve is 
a monotonous function of g reaching its maximum at g* = 1. 
Decrease of the tunnel probability leads to development of a 
well- pronounced maximum. In the large-nR limit, it is posi- 
tioned at g, — F, see Eq. (fl2l . 



Here, k € (1,til) counts the rows of an x matrix 
under the sign of determinant. Equations ([5]) and (fTU|) 
represent our second main result. 

Distribution of Landauer conductance. — Although the 
central result of this Letter, Eq. ©, allows us to address 
the problem of conductance fluctuations in full generality, 
the most explicit formulae can be obtained for the illus- 
trative example of an asymmetric cavity whose left (non- 
ideal) lead supports a single propagating mode (til = 1). 
For such a setup, a probability density /(i, nR )(ff;r) of 
the Landauer conductance is proportional to the mean 
density i^( l nR )(i?, R) of reflection eigenvalues taken at 
R = 1 — g. Materialising this observation with the help 
of Eq. (Uni), we derive^: 



/(i, I1R )(5;r) = / (1 ,„ R )(g;i)r"-+ 1 

x 2 F 1 (n R + l,n R + l;l; (1-T)(1 



9))- (11) 



Here, T is the tunnel probability of the left point contact, 

describes the conductance 



whilst /(i,„ R )(g; 1) = n R g n * 1 
density when the left point contact is ballistic. 

The probability density of Landauer conductance 
Eq. (fTTj) shows an unusually rich behavior (see Fig. [!}. 
First, it exhibits a pronounced maximum whose position 
<?*, for a generic value of the tunnel probability T, de- 
pends on the number (riR,) of propagating modes in the 
ideal lead. Second, numerical analysis of Eq. (fTTj) reveals 
existence of a "critical" value (Fo) of the tunnel proba- 
bility: for r < To, increase of 71r makes the maximum 
position move from left to right until it approaches its 
saturated location g* — T; on the contrary, for T > To, 
as tir increases, position of the maximum moves in the 
opposite direction eventually reaching g* = T. 

To describe this effect analytically, one has to seek an 
explicit functional form of g*(r,nR,) for arbitrary F and 



«r, which appears to be an impossible task. However, 
some progress can be made in the large-riR, limit, when 
an l/riR expansion can be developed. A somewhat cum- 
bersome calculation^ based on the asymptotic analysis 
of the hypergeometric function in Eq. (|11[) brings out the 
remarkable formula 



3*(r,n R ) = r 



1 



2n R V 7 



0(n 



R > 



(12) 



suggesting that the "critical" value To of the tunnel prob- 
ability equals Tq — 6/7. This prediction is unequivo- 
cally supported by numerics based on the exact Eq. (|lll) , 
see Fig. [2j We believe that experimental testing of the 
"6/7" effect may be feasible within the current limits of 
nanotcchnology-2£. 

Finally, we mention that a calculation of /( nL! „ R ) (<?; T) 
becomes increasingly complicated for > 1. This diffi- 
culty, however, can be circumvented by focussing on the 
moment / cumulant generating function™ of the Landauer 
conductance, that can be related (under certain assump- 
tions) to solutions of the two-dimensional Toda Lattice 
equation^. 

Sketch of the derivation. — Having discussed a few (out 
of potentially many) implications of the joint probability 
density of reflection eigenvalues in chaotic cavities probed 
via both ballistic and tunnel point contacts [Eq. ©], 
let us outline its derivation. The (3 = 2 Poisson ker- 
nel Eq. ([2]) with the average scattering matrix S set to22 
S = diag(7L,0 x 1„ R ) and a polar-decompose d 14 i 21 uni- 
tary scattering matrix S, 

Ml 

v x 



u 2 
v 2 



(13) 



is our starting point. Here, 



AA T 



iX 



(14) 




FIG. 2: Unnormalised conductance "distribution" 
n R 1 /(i.n R )(ffi r°) plotted for the "critical" tunnel prob- 
ability Fo = 6/7 and various numbers of propagating modes 
riR in the ideal lead: wr, = 10 (solid line), iir = 20 (long- 
dashed line), = 30 (dashed line), iir = 50 (dotted line), 
and riR = 100 (dot-dashed line). Position of the maximum is 
almost "frozen", depending very weakly on rr [see Eq. (|12[l ]. 
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the matrix A is an x zjr rectangular diagonal ma- 
trix such that AA T = diag(T\, . . . , T„ L ) if < ur,, and 
AA T = diag(Ti, . . . , T„ R ; 0x1 

m-na) otherwise; the ma- 
trices Uj and Vj are unitary matrices of the size x 
and 71r x riR, respectively. Restricting ourselves to a 
structurally more transparent case 2 ^ til < tir,, we notice 
that the polar decomposition induces the relation 



dM(s) = p (in L - h) n ^ n d ^ u ^ Mv a ), (15) 

3=1 a=l 

where Po(l nL — R) = Pq(T) is the joint probability den- 
sity function of transmission eigenvalues at j3 = 2 in case 
of ideal leads [Eq. ([6])]. and dfj, is the invariant Haar mea- 
sure on the unitary group. 

Substituting Eqs. (HSJ) and flU]) into Eq. © taken at 
/3 = 2, and considering the elementary volumes identity 
Eq. (p~5|) . we conclude that the j.p.d.f. of reflection eigen- 
values in the non-ideal case admits the representation 

P m o)(R) oc det N (t n ^-jl)P (t n ^-R) 
x / dfx(U) I dn(V) 

iu(m) iu(n L ) 

x dct- JV (l nL -TtE/ffVt) 

x det- JV (l„ L -y^C/HL), (16) 



where the notation q stands for g = \J 1„ L — AA T = 

diag(i?^ 2 , . . . , Rni 2 ). Notice, that for any finite 7l, the 
U(ul) xU(til) group integrals in Eq. (fT6|) effectively mod- 
ify the interaction between reflection eigenvalues, which 
is no longer logarithmic [see Eq. ([5])]. 

The group integrals in Eq. (|16|) can be evaluated by 
employing the technique of Schur functions 33 and the the- 



ory of hypergeometric functions of matrix argumen t 34 ' 35 . 
(An alternative derivation, based on the theory of r func- 
tions of matrix argument, was reported in Ref.—.) Leav- 
ing details of our calculation for a separate publication 2 ^ , 
we state the final result: 



det 0\fc)e(l,n L ) 2^1 («R + 1,«R + 1; 1; IjRk) 



A„ L ( 7 2 )A„ L (#) 



(17) 



Combining the last two equations together, we reproduce 
the joint probability density function of reflection eigen- 
values announced in Eq. ©. 

Summary. — In this Letter, we have outlined an RMT 
approach to the problem of universal quantum transport 
in chaotic cavities probed through both ballistic and tun- 
nel point contacts. While our central result Eq. ^ marks 
quite a progress in equipping the field with nonpertur- 
bative calculational tools, certainly more efforts are re- 
quired to bring the theory to its culminating point: (i) 
relaxing a point-contact-ballisticity for the second lead, 
(ii) extending the formalism to other Dyson-Altland- 
Zirnbauer symmetry classes! 3 -' 3 ^, and (iii) studying in- 
tegrable aspects of the theory, much in line with Ref. 31 , 
is just a partial list of related challenging problems whose 
solution is very much called for. 

This work was supported by the Israel Science Foun- 
dation through the grant No 414/08. 

Note added in proof. — Recently, we have learned about 
the paper by Y. V. Fyodorov^ who studied a probability 
of no-return in quantum chaotic and disordered systems. 
In a certain limit, this probability can be reinterpreted 
as the Landauer conductance distribution /(i,ri R )(<?; T) 
given by Eq. ((11)) of this Letter. We have explicitly veri- 
fied that both results are equivalent to each other. 
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